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ABSTRACT 


This  paper  deals  with  systems  of  singularly  perturbed  ordinary 
differential  equations  posed  as  boundary  value  problems  on  an  infinite 
intervals  The  system  is  assumed  to  consist  of  singularly  perturbed 
(fast)  components  and  unperturbed  (slow)  components  and  to  have  a 
singularity  of  the  second  kind  at  Under  the  assumption  that  there 

is  no  turning  point  we  derive  uniform  asymptotic  expansions  (is  the 
perturbation  parameter  tends  to  zero)  for  the  fast  and  slow  components 
uniformly  on  the  whole  infinite  line.  The  second  goal  of  the  paper  is 
to  derive  convergence  estimates  for  the  solutions  of  'finite'  singular 
perturbation  problems  obtained  by  cutting  the  infinite  interval  at  a 
finite  (far  out)  point  and  by  substituting  appropriate  additional 
boundary  conditions  at  the  far  end.  Using  a  suitable  choice  for  these 
boundary  conditions  the  order  of  convergence  is  shown  to  depend  only  on 
the  decay  property  of  the  infinite  solution.  4 _ 
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SINGULAR  PERTURBATION  PROBLEMS 
WITH  A  SINGULARITY  OF  THE  SECOND  KIND 


PETER  A.  MARKOWICH+  AND  CH.  A.  RINGHOFER* 

1.  INTRODUCTION 


In  this  paper  we  deal  with  the  singular 

(1.1)  Ey'  «  t°h(y,z,t,£) 

(1.2)  z*  *  tag(y,z,t,e) 

(y)  e  c( [i,»] )s  <—>  (y)  e 

2  2 


a  > 


-i,  t  e  (i,-) 


C([1,*))  and 


Nris  gra&i 
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(1.3) 


(1.4) 


linyy(t,£)%  =  cy(“#e ) \ 
^z(*,e J ' 


is  finite 


«••(«::!})  -  »<*> 


where  0  <  £  <<  1» 
component  and  z 
0(E)  e  Rk  where 

(1.5) 


y,h  are  n-vectors  z,g  are  m-vectors.  y  is 
is  called  slow  component.  F(£)  is  a  k  *  (n+m) 
k  <  n+m  holds  if  the  matrix 


_uL 

•(? 


(y(-.E)/*(-/  c)  ,«,E) 


called  fast 
matrix. 


has  at  least  one  eigenvalue  with  positive  real  part.  For  o  >  1  the  system 
(1.1),  (1.2)  has  a  singularity  of  the  second  kind  of  t  =  °». 
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For  this  we  use  techniques  already  developed  for  'finite*  singular 
perturbation  problems  as  for  example  matched  asymptotic  expansions  (see 


O’Malley  (1978),  (1979),  Ringhofer  (1980),  (1981))  and  the  theory  of  singular 
boundary  value  problems  (see  de  Hoog  and  Weiss  ( 1980a, b),  Markowich 
(1980a, b,c)  and  Lentini  and  Keller  (1980)). 

We  show  that  the  solutions  y,z  of  (1.1),  (1.2),  (1.3),  (1.4)  fulfill 

(1.9)  z(t, e)  =  z(  t )  +  0(E),  t  e  [1,“] 

(1.10)  y(t,e)  =  )  +  y(t )  +  0(e),  t  e  [1,-] 

where  y,z  are  solutions  of  (1.1),  (1.2),  (1.3)  with  e  «  0  (reduced 
problem)  fulfilling  appropriate  boundary  conditions.  Here  ®(T)  decays 
exponentially  to  zero  as  t  ♦  "  (boundary  layer  term)  and  z,y  decay  to  a 
finite  limit  as  t  +  *•  fulfilling 

(1.11)  0  =  h(y-,zoo,“,0) 

d.i2)  o  =  g(y„,z„,*#0). 

This  result  generalizes  the  results  by  O'Malley  (1979)  and  Ringhofer  (1980) , 
(1981)  obtained  for  finite  interval  singular  perturbation  problems. 

Singularly  perturbed  initial  value  problems  on  the  infinite  line  have 
been  investigated  by  Hoppensteadt  (1966)  by  imposing  severe  stability 
assumption  on  the  reduced  problem. 

The  second  goal  is  to  study  approximating  'finite*  singular  perturbation 
problems,  which  are  set  up  by  cutting  the  infinite  interval  (1,*]  at  a 
finite  point  T  »  1  and  by  substituting  (for  the  continuity  condition  (1.3) 
at  t  =  “)  additional,  so  called  asymptotic  boundary  conditions  obtaining  a 
'finite*  singular  perturbation  problem 


(1.13) 


eyj.  -  t  h(yT,zT,t,e) 


1  <  t  4  T 


(1.14) 


^  *  t  g(yT.zTt,e) 
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a  i  aii  i^ 


(1.15) 


(1.16) 


F(C  )( 


yT(1,e,1 
*T<1 'C) 


8(e) 


s(t,e  )( 


yT(T.£) 

VT'e) 


Y(T,e) 


where  S(T,e)  is  an  (n+m-k)  x  (n+m)  matrix,  Y(T,e)  e  Rn+m  The  condition 

(1.16)  shall  reflect  the  asymptotic  behaviour  of  the  'infinite'  solution 
(y,z)  as  t  ♦  “. 

'Finite'  approximating  two  point  boundary  value  problems  (for  unperturbed 
infinite  problems)  have  been  studied  extensively  by  de  Hoog  and  Weiss  (1980a), 
Markowich  (1980b)  and  Lentini  and  Keller  (1980). 

We  show  that  under  rather  mild  assumptions  on  the  'infinite'  problem  (a 
certain  'wellposedness'  is  required)  there  is  a  choice  of  S(T,c)  =  s  and 
Y(T,£)  =  Y  only  depending  on  the  reduced  (£  =  0)  infinite  problem  such 
that  the  'finite'  (perturbed)  problem  has  a  unique  solution  yT,zT  for  T 
sufficiently  large  and  £  sufficiently  small  (but  T  and  e  independent) 
which  fulfills  the  convergence  estimate 

(1.17)  ^z-zT^[1,T]  <  const  ( exp( -  T°+1)  +  e) ,  c  >  0 

(1*1 .  denotes  the  sup-norm  on  [a,b] )  where  the  constants  are  independent 

la  ,b] 

of  T  and  £  . 

The  'finite'  singular  perturbation  problem  (1.13),  (1.14),  (1.15),  (1.16) 
can  then  be  solved  by  polynomial  collocation  methods  (see  Kreiss  and  Nichols 
(1975),  Ringhofer  (1981),  Ascher  and  Weiss  (1981)).  An  exponential  mesh  size 
strategy  for  'long  interval'  problems  has  been  developed  for  the  Box-scheme  by 
Markowich  and  Ringhofer  (1981).  This  can  be  used  on  [<*»,T] ,  w  >  1  while 
within  the  boundary  layer  (on  [1,1+0(e |ine| ) J)  a  very  fine  grid  (see  Ascher  and 


Weiss  (1981))  has  to  be  used.  Since  the  solution  of  (1.13),  (1.14),  (1.15), 
(1.16)  is  smooth  (has  uniformly  (in  C)  bounded  derivatives)  on  [1+0(e]£ne|) ,oj] 
standard  techniques  can  be  used  there. 

The  paper  is  organized  as  follows.  In  chapter  two  linear  constant 
coefficient  problems  are  treated,  in  chapter  three  variable  coefficients  are 
admitted  and  chapter  four  is  concerned  with  nonlinear  problems. 
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2.  Constant-Coefficient  Problems 


At  first  we  study  the  problem 


(2.1) 


€y*  =  taA(c)y  +  t°B(e)z  +  taf(t,£)''' 


(2.2) 


=  taC(e)y  +  taD(e)z  +  taq(t,£) 


1  <  t  <  « 
a  >  -1 


(2.3) 


F(e  )( 


y(1*e  )■> 
z(1*c  V 


6(c ) 


(2.4) 

(;)  «  C(  [1,®] ). 

Here  y 

,f  are  n-vectors,  z  and  g 

are 

m-vectors,  A(e) 

is 

an  nxn- 

matrix. 

B(e )  an  n^m-matrix,  C(e ) 

an 

mxn-matrix  and 

D(e  ) 

is  an  fflxm« 

matrix. 

The  dimension  of  the  matrix 

F(e) 

and  the  vector 

0  (£ ) 

will  be 

discussed  in  the  sequel  (obviously  F(e )  has  n+m  columns).  We  assume  that 
(2. .5)  f,g  6  C((1/»]  x  (Ofeol),  eQ  >  0 

(2.6)  A,B,C,D,F,0  G  C([0,eol). 

and  that  A,B,C,D, f ,q,F,0  are  (uniformly)  Lipschitz  continuous  at  e  *  0. 

It  is  convenient  to  decouple  the  system  (2.1),  (2.2)  by  a  linear  trans¬ 
formation  such  that  the  unperturbed  equation  of  the  transformed  system  does 
not  contain  the  fast  component.  We  use  the  transformation  given  by  Ascher  and 
Weiss  (1981): 


<2-7’  ft?- 

Assuming  that  A  1  (e )  exists  for  e  G  tO.Ep)  and  that  lA_1(e)l  <  const, 
for  £  G  (0,£g]  we  determine  L(c )  from  the  equation 

L(e)  -cfcjA^fe)  -  e(D(e)L(e)A'’1  (e)  +  L(e )B(e )L(e )A_1 (e ) ). 


(2.8) 


A  simple  contraction  argument  assures  the  unique  solvability  of  (2.8)  such 


that 

(2.9) 

L(e)  =  c(e)A-1(e)  +  0(e) 

holds. 

The  new  system  (with  u,v  as  dependent  variables) 

(2.10) 

eu*  =  taA(6)u  +  t°3(e)v  + 

t°f (t,e) 

(2.11) 

v*  =  taD(6)v  + 

t°g(t,ej 

(2.12) 

(“)  e  C( [1,«] ) 

where 

A  =  A  +  EBL,  5  =  D  -  LB,  g  =  g  -  Lf 

holds. 

Now  we  make  assumptions  on  the  eigenvalues  of  A(0),  0(0)  **  D(0)  - 

-C(0)A~1(0)B(0): 

A( 0 )  has  r+  eigenvalues  with  positive  real  part  and  r_ 
(2.13)  eigenvalues  with  negative  real  part  (counting  algebraic 

multiplicities)  and  r+  +  r_  =  n. 


(2.14)  5(0)  has  r+  eigenvalues  with  positive  real  part  and 

r_  eigenvalues  with  negative  real  part  and  r+  +  r_  =  m. 
We  investigate  the  perturbed  problem 


(2.15) 

eu; 

*  taA(0)u1  +  t°B(e)v1  +  t^UfC) 

a-  o 

(2.16) 

v; 

(U’' 
'■v  ; 

*  t  D(0)v1  +  t  g^(t#c) 

(2.17) 

1  e  c( [i,-i ). 

V1 

The  assumptions  (2.13),  (2.14)  guarantee  that  there  are  transformation 
matrices  Ej,  E2  such  that  the  matrices  Jj,  J2  defined  by 
(2.18)  (a)  A( 0 )  -  E^e"1  (b)  5(0)  -  E^e"1 

are  in  block  diagonal  form: 


(2.19) 


,  f*  _+  _  ^ 

JO  J2  0  ! 

(a)  J.  -  ,  (b)  J  * 

o  j;  0  J 


r+ 


where  only  have  eigenvalues  with  positive  (negative)  real 

parts.  We  substitute 


(2.20) 


u„  E.  0 

(.)-C  e2HP 


and  get 

a  or- 

(2.21) 

£w‘  *  t  +  t  B  ( E )  x  + 

a 

(2.22) 

X*  »  t  JjX  + 

(2.23) 

(*)  e  C(  {1,-D 

where  B(E)  -  e"1B(C)E2,  -  e"1^,  *  e"^  holds.  We  solve  (2.22), 

(2.23)  defining 


(2.24) 


a+l_£<x+1 

♦  (t,fi)  -  exp(” — -  J2) 


such  that  for  the  slow  component 


(2.25) 


x(t,E)  -  $(t,1)  w  £  +  (H^C.elKt) 
r 


holds  where  (  e  C  *  and  the  solution  operator  for  6  >  1  is  defined  by 

...  .  .  1  .  .  I  _  ft  .  ,  .  f  .  _  Xi.®  /.  r*  \  J—  x 


(H5g(»,E) )(t)  -  /t^(t,5)D+#~1(s,5)s  g(s,E)ds  + 


(2.26) 


+  /*♦( t,$)D  $  1(s,6)s°g(s,e)ds. 

6 
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Here  D+,D_  are  diagonal  projections 


(2.27) 


(a)  D 


+  L 


,  (b)  D  = 


Hg  was  used  by  de  Hoog  and  Weiss  ( 1980a, b)  and  they  showed 

(2.28 )  (a)  :  C<  £« ,“] »  ♦  C( {«,•]>,  5  >  1 

(2.28) (b)  [5  »]  *  const,  independently  of  6 


(2.28) (c) 


(H.f (•,£))(»)  «  -j  f(“,£). 


Estimates  of  the  asymptotic  behaviour  of  (H^f)(t)  as  t  *  °*  are  given  in 
Markowich  ( 1980a ) . 

Inserting  (2.25)  into  (2.21)  we  can  regard  h(t,£)  “  B(£)x(t,£)  + 


+  f  ^ (t , £  )  as  inhomogenity  and  are  left  with  solving  for  the  fast  component: 


(2.29) 


£w'  *  t°J  w  +  t°h(t,£),  w  e  C( [1,®] ). 


We  define 


(2.30) 


o+l  .a+1 

-Ht,5,£)  =  exp(  J,) 


(2.31) 


where 


(Gr  *hC,£))(t)  -  /t^(t,5,£)D>"1(s,6,e)s0th(s,£)ds  + 
00  * 


+  /%< t,«,£)D  ♦-1(s,6,£)s<Xh(s,£)ds 


(2.32) 


(a)  D, 


(b)  D  = 


~  0 

I 

r_ 


holds.  Then  the  solution  of  (2.29)  is 


(2.33) 


w(t,e)  -  <|>(t. 


1,0  \  C  +  (GEf1hC,£))(t| 

-  r- J 


for  ?  e  c  . 
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It  is  easy  to  show  that  G£  ^  has  the  following  properties: 
(2.34)  (a)  Ge^  :  C(  («,-))  ♦  C(  «  >  1 

(2.34}(b)  lGe  ^1  ^  „ j  *  const,  independently  of  5,e 

(2.34)  (c)  (G£  5hC,e))(*)  -  -J^hC^E). 


If  h( *,e ),  ---af£-  ec([1,«])  for  all  e  sufficiently  small 


then 


( 2.34 )  (d) 

where 

(2.35) 


(Ge,«h<*'E,)<t)  *  -J^Wtfe)  +  *<t,«/e)D_jj1h(«,e) 
+  E(Ae^h(*,e))(t) 


“s/C'*!  \s  -]  <  const(lh(*,e)l  .  ♦  max 

'  se[5,-]  s 

holds.  In  the  sequel  we  use  the  space 


(2.36)  [«,-])  -  C(  [«,“])  n  c\ [«,-))  n  {f|  max 


se[6,»)  s 


and  as  norm  we  take  Ifl®.  s-  lf»r.  +  max  lf,(s)/sa«, 

[6'  1  IC'  5  se[6  -] 


We  rewrite  (2.11)  obtaining 
(2.37) 

and  get  by  integration  (using  (2.20),  (2.25)) 


v’  -  t°D(0)v  +  ta(D(E)  -  D(0))v  +  t°g(t,E),  t  >  1 


v(t,E)  -  E2*<t,1) 


(2.38) 


I~ 

r. 


■1 


5  +  ^2(8^2  (D(E)  -  D(0  )  )v(  *,E) )  (t)  + 


+  E2(H1E21gC,e))(t). 


where  H^sCdl,®])  +  Ca(  [!,»])  is  bounded. 
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We  conclude  from  (2.28)  (b)  that  (1  -  E^H^E^fDie)  -  5(0) )-1  =  I  +  0(e)  as  an 


operator  on  C([1,°°])  such  that 

I  0 


v(t,e  )  =  E2<J>  < t ,  1  )  ^ 


(2.39) 


Lr-. 


C  +  0(e)Ei  +  E2(H1E21q(*  ,e  ))(t)  + 


+  0(el5(*,e)l{i  wl) 


holds  uniformly  on  (1,“].  (The  asymptotics  for  v*(t,e)/t  follow 
immediately).  Rewriting  (2.10)  gives 

(2.40)  eu'  =  taA(0)u  +  ta(A(e )  -  a( 0)  )u  +  taB(e)v  +  taf(t,e),  t  >  1 


and  we  obtain  (using  (2.20),  (2.33)) 


u(t,e  )  =  E1'|>(t,  1  ,e  ) 


(2.41) 


■  on 

I 

r 

i_  si 


c  +  E1(Ge  1B“'1  (A(e  )  -  A( 0)  )u(*  ,e )  (t)  + 


-1 


+  E1(Ge  1 Ej  (B(e )v( * ,e )  +  f(*,e)))(t). 


From  (2. 34) (b)  we  conclude  that  (1  -  EG  E*1 (A(e )  -  A(0)))-1  *  I  +  0(e)  as 

1  C  f  1  I 

operator  on  C([1,*])  such  that 

(2.42)  u(t,e )  .  E1<»(t,1,e)^i0Jc  +  0(e)c  + 


+  E1(Ge(1E;1(B|t)v(.,e)  +  f  (•  ,e  ) ) )  (t)  + 


+  o(e(lvC,e)»n  +  ,f(*'c)l[i,-])* 

Resubstitution  in  (2.7)  qives  y,z.  Since  y,z  depend  on  r  +  r 
parameters  (5  and  C)  we  assume  that  the  matrix  F(e )  as  in  (2.3)  has 
r  +  r  rows.  By  collecting  the  terms  of  y,z  which  depend  on  ? ,C  we  get 
Theorem  2.1.  Let  f,g  G  cj((1,*°])  uniformly  for  small  e  >  0  and  assume 


that  the  (r  +  r  )  *  (r  +  r  ) -matrix 


(U 


-1-1  C  0 

E1D+E1  A  (0)E2  I~ 
r 


(2.43) 


is  nonsingular.  Then,  under  the  given  assumption  on  A,  B,  C,  D,  F, 8  the 
boundary  value  problem  (2.11,  (2.2),  (2.3),  (2.4)  has  for  all  e  sufficiently 


r_+  r_ 

small  and  for  all  8(£)  €  R  a  unique  solution  (^)  which  depends 

uniformly  continuous  (in  £)  on  5(e)  and  on  f,g  e  C( [ 1 ,°° ] )  when  regarded 
as  dwelling  in  C([1,*J).  Moreover  (^)  e  C^([6,“]l  for  6  >  1  depends 
uniformly  continuous  on  f,g  e  c\(1,«*J)  and 


(2.44) 


z(t,€)  -  E2<Mt#  1 ) 


[■y  • 


!2<H1E21(gl*'0) 


-  C(0)A  (0|fC,0|)|(t|  +  0(E) 


(2.45)  y(t,e  )  *  E^t 


B(e)z(i,e)  + 


+  E~^f(1,€)))  -  A~\o)(B(0)z(t,0)  +  f(t,0))  +  0(E) 


holds  uniformly  on  [1,**]. 


From  (2.28) (c),  (2.34)(c)  we  derive 


(2.46)  z(~,e)  -  —  ( D  <  0 )  -  C(0)A-1(0)B(0)  )_1(g(-,0)  -  C(0  )A-1  (0  )f  («»,0  ) )  +  0(E) 


y(«,E)  -  -A  (0 ) (B(0 )z(",0 )  +  f(-,0|)  +  0(E). 


(2.47) 


The  first  term  in  (2.45)  is  the  exponentially  decreasing  boundary  layer 
contribution  (the  thickness  of  the  boundary  layer  is  0(e|fcne|)  and  the  second 
term  is  the  solution  of  the  reduced  problem  (2.1)  (with  e  =  0). 

If  we  drop  the  restriction  that  D(0)  =  D(0)  -  C(0)A  ^(0)B(0)  has  no 

eigenvalue  on  the  imaginary  axis  (see  (2.14))  we  have  to  assume  that 

f (t,c),g(t,e)  converge  to  zero  algebraically  as  t  ♦  A  sufficient  order 
-(0+1 )r-Y 

of  decay  is  t  ,  where  r  is  the  dimension  of  the  largest  Jordan 

block  of  0(0)  which  has  an  eigenvalue  on  the  imaginary  axis  and  Y  >  0  (see 
Markowich  (1980a)).  For  the  contraction  arguments  algebraically  weighted 
C([1,-])  resp.  cj([1,-])  spaces  have  to  be  used. 

For  the  numerical  solution  of  (2.1),  (2.2),  (2.3),  (2.4)  we  cut  the 
infinite  interval  at  a  finite  point  T  >>  1  and  replace  the  continuity 
requirement  (2.4)  by  r+  +  r+  boundary  condition  at  t  =  T.  These  boundary 
conditions  shall  reflect  the  asymptotic  behaviour  of  y,z  as  t  ♦  00 .  So  we 
get  the  'finite'  singular  perturbation  problem 


(2.48) 


£y;  -  taA<e)yT  +  taB(e)zT  +  t°f (t,e) 


1  <  t  <  T 


(2.49) 


z'T  =■  tac(E)yT  +  t°D(e)zT  +  tag(t,e> 


(2.50) 


(2.51) 


rVT'eS 

S(T'E,(zT<T,e,) 


~  r++r+ 

Here  S(T,e)  is  an  (r+  +  r+)  *  (n+m)-matrix  and  Y(T,e)  e  R 


A  possible  choice  is 


(2.52) 


(2.53) 


Y(T)  5  Y(T,e)  = 


[lr  ,0]E"1A"1(O)f(T,O) 


-  [I?+']E21z(0 


,0) 


This  'asymptotic '  boundary  condition  has  been  used  by  de  Hoog  and  Weiss 
(1980a),  Markowich  (1980b)  and  Lentini  and  Keller  (1980)  for  unperturbed 
problems  on  infinite  intervals.  Let  J^t,e)  denote  the  fundamental  matrix  of 

(2.1) ,  (2.2)  (Jt(l,e)»  I),  Then  by  proceeding  as  de  Hoog  and  Weiss(1980a)  did 

we  can  easily  show  that  sMt,6)  does  not  contain  exponentially  decreasing 
terms.  Therefore  the  countary  condition  *  0  sets  the  expo¬ 

nentially  increasing  solution  components  of  the  homogenous  problem  (2.1), 

(2.2)  to  zero.  Y(T)  as  in  (2.53)  is  the  necessary  (boundary)  correction  term 
for  the  inhomogenous  problem. 

By  proceeding  as  in  de  Hoog  and  Weiss  (1980a)  we  find  the  stability 
estimate  for  the  solution  of  (2.48),  (2.49),  (2.50),  (2.51)  when  using  S  as 
in  (2.52)  and  assuming  that  the  matrix  (2.43)  is  nonsingular 


(2.54) 


<  const ( «  8(e ) I  +  lY(T,e)l  +  »fC,e)l  + 

L  * 


+ 


lgC,e)l 


[1,T] 


) 
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r  +r  r  +r 

for  all  0  6  R  ”  ~,  y  6  R  +  +  ,  f,g  G  C([1,T])  where  the  constant  is 
independent  of  T  and  e.  By  subtracting  (2*48),  (2.49),  (2.50),  (2.51) 
from  (2.1),  (2.2),  (2.3)  we  get  the  error  estimate 

<2-55> 

Inserting  (2.44),  (2.45)  into  the  right  hand  side  of  (2.55)  (when  using 
(2.52),  (2.53))  gives  the  uniform  estimate 

y-y 

if  „*)■»,  «,  <  const(lE  (H  E~* (g (•  ,0)  -  C(0)A\o)f  (•  ,0) )  (t) 
(2.56)  T 

-  *(-,0)1  +  0(e)). 

Convergence  follows  because  (2.28 )(c). 

Since  the  solutions  of  the  reduced  problem  (e  =  o)  (2.48),  (2.49)  do 
not  generally  fulfill  (2.51)  one  has  to  expect  a  boundary  layer  at  t  «  T 
whose  height  can  be  estimated  by  the  right  hand  side  of  (2.56). 

Estimates  of  the  order  of  convergence  of  the  first  term  of  the  right  hand 
side  of  (2.56)  depending  on  the  decay  of  f,g  as  t  ♦  *»  are  given  in 
Markowich  (1980a),  (1980b). 

Therefore  under  the  given  assumptions  the  asymptotic  boundary  condition 
(2.51)  can  be  consructed  with  respect  to  the  reduced  (e  =  0)  infinite 
problem. 
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3.  Variable  Coefficient  Problems 


We  consider  the  problem 
(3.1)  €y* 


t°A(t,e)y  +  t°B(t,e)z  +  t°f(t,e)^ 


> 

(3.2)  z*  *  t°C(t,e)y  +  t°D(t,e)z  +  t°g(t,e) 


1  <  t  <  * 
a  >  -1 


(3.3) 


Fiei  ■  »<*> 


(3.4)  (y)  e  C(  [1,-n 

Z 

where  the  dimensions  are  as  in  Chapter  1  and  assume  that 

(3.5)  A,B,C,D,£,g  e  C( [1,-J  x  [0,6^,  F,S  «  C( [0,eQ) ) 

holds  for  some  >  0  and  that  F, fl,A,B,C,D,f ,g  are  uniformly  Lipschitz 
continuous  at  e  =  0. 

Moreover  we  assume  that  the  eigenvalues  X(t )  of  A(t,0 )  split  up  into 
two  groups  such  that 

(3.6)  Re  X^t)  >  c+,***,Re  Xy  (t )  >  c+,  c+  >  0,  t  >  1 


(3.7)  Re  X  (t)  <  -c  ,**,,ReX  (t)  <  -c  ,  c  >  0,  t  >  1  (n  -  r  =»r) 

r++i  -  n  -  -  +  - 

holds  (eigenvalues  are  counted  accoring  to  algebraic  multiplicities)  and  that 
there  is  a  transformation  to  block  form 


A(t# 0)  -  E(t)J(t)E-1(t),  J(t) 

such  that  the  eigenvalues  of  J+(t) (J_(t) 

(*  ,.(t ),••*, A  (*)>  and 

r.Ti  n 


J  (t) 

0 

+ 

0 

J_(t) 

_ 

_ 

are  Aj(t),«##,A  (t) 


(3.9) 


I  El 


Mf-] 


<  const . 
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Under  the  assumptions  (3.6),  (3.7)  and  additional  smoothness  assumption  on 
A  this  transformation  matrix  E  exists  at  least  locally  (everywhere  in 
[1/*]),  but  the  assumption  of  the  global  existence  is  much  more  restrictive 
(see  O'Malley  (1979)).  At  first  we  investigate 

(3.10)  ey'  -  t°A(t,e)y  +  y  e  C({1,»)) 

and  substitute 

(3.11 )  y  *  E(t)x 
obtaining 

(3.12)  ex'  -  t°J(t)x  +  ta(E-1(t)(A(t,e)  -  A(t,0))B(t)  - 

-et"aE"1(t)E'(t))x  +  taE-1(t)h(t,e) 

(3.13)  x  e  C( [1,“]). 

Using  a  perturbation  approach  we  first  solve 

(3.14)  £u'  -  taj(t)u  +  t0d(t,e),  u  e  C([1,°*]). 

According  to  (3.8)  the  system  (3.14)  splits  up  into 

(3. 14 )  (a)  eu^  »  t“j+(t)u++  t°d+(t,e),  u+ec([1,-]) 

(3.141(b)  cu'  -  t°J_(t)u_  +  t°d_(t,£),  u_  e  C([1,»]). 

At  first  we  analyse  (3.14a),  which  we  rewrite  as 

( 3. 15 )  (a )  eu|  -  taj+(-)u+  +  t°(J+(t)  -  J+(»))u+  +  t°d+(t,c) 

( 3. 15 )  (b )  u+  e  C< [1,*] ). 

We  regard  (3.15)  as  an  inhomogeneous  constant  coefficient  problem  with  the 
fundamental  matrix 

(3.16)  *+<t,«,e)  -  •*p(^7,(ta+1-  «a+1)),  «  >  1 

and  with  solution  operator 

(3.17)  (G^6d+C,e))(t)  -  ^  /t*+(t,«,cH»*1(8,4,e)sad+(s,€)ds,  t  >  6. 
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•ftie  solution  of  (3.15)  is 


(3.10)  u+  =  G^6(J+C)  -  J+(-))u+  +  Ge  ,6d+(*  ,e  >* 

Since  (2.39)  holds  and  J+(t)  ♦  J+(*>  the  operator  I  -  Ge  fi(J+(*)  -  J+(")) 
is  invertible  on  C<  (6 ) )  for  5  sufficiently  large.  We  obtain 

(3.19)  u  (t)  =  ((I  -  6(J+C)  -  J+(-)))_1Gg  6d+C,e))(t),  t  >  «. 

^ - r— - ' 

(9*d  (• »e ) ) (t) 

c  + 

Tto  qet  a  solution  on  (I.*!  we  solve  the  termif*^  problem 

(3.20)  eu|  *  t°J+(t)u+  +  t°d+(t,e),  "*  '*  ■  '  '3 

(3.21 )  u+(«  )  -  u+(6  ). 

and  set  (0+d  (•  ,e  )  )<t)  :=  u  (t)  for  t  <  I.  *  is  an  operator  on  (1,-1 
€  +  ▼  c 

and  since  the  eigenvalues  of  J+(t)  have  strictly  positive  real  part 

(3.22) 

holds  and  because  of  (2. 34d),  (2.35) 

i«0 


,0e'[  1,-]  «  ConSt 


(9+d+e,e))(t)  =1  ( (J+(* )  -  J+(-)))V6d+C,e))(t) 


(3.23) 


-  -{J+(t))_1d+(t,c)  +  0(eld  +  e,e)l~6<1D]) 
holds  for  t  >  6.  By  continuation  (3.23)  holds  for  t  >  1  (See  Ringhofer 


(1981)). 


We  rewrite (3. 1 4) (b)  analoguously 

( 3. 24)  (a)  eu‘  =  t\j_(«)u_  +  t“{J_(t)  -  J_(-))u_  +  t°d_(t,e) 

( 3. 24)  (b)  u  6C((1,-1) 


and  define  the  fundamental  matrix 

.  ^  ,^a+1  ,a+1,| 

(3.25)  »_(t,6,e)  =  exP^-i7STTP (t  "  6  )>r 


6  >  1 


and  solution  operator 

(3.26)  (G~,3  (•»£)(  t)  ■  —  J  (t,6,e)»|>  1  (s  ,6  ,e )  j  d_(s  ,c  )ds, 

e,s  -  e  g  ~ 

such  that  the  general  solution  of  (3.24)  is 


t  >  6 


(3.27) 


u_  =  (I  -  G“  5(J _(•)  -  J_(»)))_14»_(»,6,e)C  + 

+  (1  -  G‘  (J_(* )  -  J_(»)))-1G£  gd_(*  ,£ ) 
r 

for  C  6  C  and  t  >  6.  We  call  the  first  term  on  the  right  hand  side  of 

(3.27)  ijJ  (t,6,e)C  and  the  second  u  (t,£).  Obviously 

P_ 

(3.28)  (6,6,e)  =  I  ,  (i  (•,«,£)  =  0,  u  (6  ,e )  =  0 

r_  p_ 

hold,  if  has  a  boundary  layer  at  5.  Ihe  homogenous  problem  (3. 14)(b)  has 
a  fundamental  matrix  ij»  (t,£)  such  that 

— c 

(3.29)  (a)  *_<1,£)  *  Ir  ,  (b)  M>_(t,£)l  <  Clexp(— (t°+1-1)) 

holds  for  t  G  [1,6]  where  c^Cj  >  0  (see  Ringhofer  (1981)  and  under  more 
general  assumptions  O’Malley  (1978)).  We  set 

(3.30)  *Jt,E)  «  *_(t,6  ,£)♦_(«,£),  t>1. 

A  * 

Since  i|>  (6  ,e  )  =  ip  (5  ,e )  we  obtain  ^  and  the  boundary  layer  has  been 

shifted  from  6  to  1.  On  [1+  0(£  |ine  | )  ,«•  ]  the  matrix  i|»  is  smooth. 
Another  particular  solution  is 

(3.31)  u  ( t  ,£ )  (t,E  )<i-1  (S,E  )s°d  (s,£)ds,  tG  [1,6], 

p  E  i 


Since  •♦_(t,£  )^_1  (s,e  )l  <  c3  exp(e^q^  ^  (ta+1-sa+1 ))  holds  on  [1,6]  we 


derive 

(3.32) 

<  constl d  (•  ,£)!.,  - 

+  [1,6]. 

Setting 

(3.33) 

A 

u  (t,5)  =  (8  d 
P_  £ 

(«,e))(t)  ♦  ( t  ,6  ,£  )u  (6  ,£  )  +u  (t,£) 

P_  P_ 

we  obtain 

A 

u  =  u  and 

P  P 

(3.34)  *®e*[1  •]  *  c°ost 

because  on  [1,6]  we  use  (3.32)  and  on  [6,«]  we  use  estimate  (3.33)  and 
(3.27).  As  general  solution  of  (3.14)(b)  we  take 
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(3.35) 


and  we  find 


(3.36) 


u_(t,e)  -  5»_(t,e)C  +  (0”d_(*,e)  )(t),  t  >  1 

u  (t,e)  -  4»_(t,e)(5  +  j"1  ( 1  )d_(  1  ,e) )  -  j”1 <t)d_(t,e)  + 
+  0(eld_C,e)l®1  ) 


uniformly  on  [I,*]. 


Setting  9  write  the  solution  of  (3.12),  (3.13)  as 


(3.37) 


-  °  1C  +  0C<®”1  <A( -,e>  -  AC, 

♦_<*#e)J 


-1. 


0)  )E  -  EB)X  +  0£E  hC,E) 


where  E(t)  -  t^E-1  (t)E' (t )  has  been  set.  (3.5),  (3.9)  guarantee  that 
A(t,€)  -  E‘\t)(A(t,e)  -  A(t,0))E  -  £E(t)  ♦  0  as  e  ♦  0  uniformly  on 
[1,-].  Therefore  (I  -  QgM* ,e) )_1  exists  on  C([1,»l)  and  is  bounded 
uniformly  in  t  such  that 

x(t,e)  -  ((I  -  (°  eJ^(t>C  + 

(3.38) 


holds  for 
(3.39) 


C  e  C  . 
y(t,e)  - 


+  ((I  -  0eAC 


v  ,ro* 

I(t,e>. 

L  r-_ 


p  -  A 


-1 


,e))"10EE”1hC,e))(t),  t  > 
(t,0)h(t,e)  +  0(«lh(«,e)l®1j-J 


1 


J. 


t  >  1 


holds  where 


is  the  boundary  layer  term  (at  t 


1)  fulfilling  the 


r 

estimate  (3.29) (b)  and  p«C  ”.  This  is  proven  by  using  the  series 
expansion  of  (3.38). 

Returning  to  the  coupled  problem  (3.1),  (3.2),  (3.4)  we  assume  that 
(3.40)  A,B,C,D,f,g«cJ(|1,»)l  uniformly  in  e. 
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From  (3.39)  we  get  for  fixed  z  e  Ca([1,*]) 


y(t,e)  »  l(t,€) 


(3.41) 


(t,0 ) (B(t,0 )z(t,e )  +  f(t,0|)  + 


+  e(L^1,z)(t)  +  e(Iig21  f )  (t) 


where  -  C^<  H  )  *  C( £1 2 )  and  lL*1,l®  „  <  const,  i  >  1,2. 

e  u  e  l 1 9  J 

Inserting  (3.41)  into  (3.2)  gives 

z‘  =  t“(D(t,0)  -  C(t,0)A-1(t,0)B(t,0))z  +  t°E(Lg3,z)(t)  + 


(3.42) 


+  ta(l44,f)(t) 


+  t°(C(t,e)X(t,e)^°  Jp  +  g(t,e ) ) , 


z  e  c( [i,«] ), 


Again  the  operators  s  C*([1,“])  ♦  C([1,“]),  j  -  3,4  are  uniformly  (in 


-1, 


e)  bounded. 

Setting  D(t,e)  -  D(t,E)  -  C(t,€)A  '(t,e)B(t,E)  and  assuming  that 
(2.14),  (2.18b)  holds  for  O(w,0)  we  can  solve 
(3.43) 


z*  =*  t°D(t,0)z  +  t°g(t,e),  z  e  C*([1,«»]> 


by  using  the  theory  developed  by  de  Hoog  and  Weiss  (1980a,b).  We  obtain  for 


€  e  c 


(3.44) 


1  I  0  i 

z  *  (IHSjH^e”' (D( *,0 )  -  5<-,0))E2)"  ♦(•,«)  5  + 

r 


+  (I-E2Hse‘1(DC,0)  -  6(-,0||B2)‘,E2H4eJ15c,I|i  t  >  « 
where  Hj,$(t,$)  are  defined  in  (2.26),  (2.24)  reap,  and  6  is  sufficiently 
large.  The  right  hand  side  of(3.44)  can  be  continued  to  [1,4]  and  we  obtain 


(3.45) 


z(t,e)  -  *<t)L  |C  +  <Ig(«,e))(t),  t  >  1 
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4 


<  const.  Applying  this  to 


where  T  ;  C<[1,*1)  ♦  <£([1,«1)  and 


(3.42)  gives 


(3.46) 


z(t  ,e)  *  $(t) 


R]" 


e(rL^3,z)(t)  +  (ri<4,f)(t)  + 


+  <rg(*,e))(t)  +  <r(CC,e)I(*,e|^J  "j)(t)p 


rLi3)  !  ♦  C^(  £1,-1 )  and  •ri/3)l“  ,,  <  const.  Therefore 

E  Gt  Gt 

(I  -  erL<3>,  1  exist  as  operators  on  C^([1,*])  for  e  sufficiently  small 


and 


(3.47) 


z(t,£)  -  4(1: 1^°  1C  +  (r(9<*'0)  “  C(S0)A*1C,0)f(»f0)))(t)  + 


+  <rc(*#«ilc#«)j^0  Jxtjp  + 

Using  the  exponential  decay  of  | 

LrJ 


0(e). 


and  the  definition  of  Hg  it  is 


easy  to  show  that 


(3.48) 


irC( 


0(e). 


So  we  obtain 

Theorem  3.1.  Let  the  given  assumptions  on  A,B,C,D,f ,g  hold  and  assume  that 
the  (r  *  r  )  *  (r  x  r  ) -matrix 


F(0) 


2(1) 


[■I] 


E(1)D+E(1)”1a“\o)B(0)^(1)["1° 


♦o>  jZ 


in  nonsingular  (F(C)  is  a  (r_  +  r_)  *  (n+m) -matrix) .  Then  the  boundary 
value  problem  (3.1),  (3.2),  (3.3),  (3.4)  has  for  sufficiently  small  e  and 


for  all  3(E)  e  R  ,  f,g  e  C^([1/®])  (uniformly  in  E)  a  unique 

solution  y,z.  The  continuity  statements  of  Theorem  2.1  hold  and 

— 

*  .  - 

z(t,E)  =  <Mt)  J  5  +  (r(gC,0)  -  CC.OIA  C,0lf(*,0))l(t)  + 

(3.49)  _  r- 


+  0(e),  5  e  c 


(3.50) 


(t,£)  =  I(t*E)|~i0Jp  -  A-1(t,0)B(t,0)z(t,0)  +  f(t,0))  + 


+  0(e),  y  e  c 

hold  uniformly  in  [I,00].  y("»e)/  z(°°,£)  are  as  in  (2.46),  (2.47)  when 

A(0),  B( 0 ) ,  C(0),  D( 0 )  are  substituted  by  A(»,0),  B(“,0),  C(°»,0),  D(~,0). 


The  'finite'  problem  is 


(3.51) 


(3.52) 


=  taA(t,£)yT  +  taB(t,£|zT  t“f(t,E) 


z^,  =  ta0(t,£)yT  +  taD(t,E)zT  +  tag(t,E) 


1  <  t  <  T 


yTo-e>. 


(3‘53)  F(E,(z  ( 1 ,  E  ) '  *  6<€> 


(3.54) 


,VT'e,1 

S(T'e’(zT(T,E)^  "  Y<T'e) 


r++r+ 

where  Y(T,£)  e  R  ,  S(t,£)  is  a  (r++r_)  x  (n+m)-matrix.  We  assume  that 
(2.18)(a)  holds  for  A(®,0).  Then  by  proceeding  as  de  Hoog  and  Weiss  (1980a) 
did,  we  find  that  we  can  use  (2.52),  (2.53)  in  order  to  set  up  the  asymptotic 
boundary  condition  (3.54).  (We  do  not  have  to  know  E(t)  explicitely  since 


it  can  be  chosen  such  that  E(t)  +  E  as  t  +  *)•  The  convergence  estimate 
(2.56)  holds  if  we  add  0(exp(  -  ■  T°+1)),  P  >  0,  which  is  an  estimate 

A 

for  the  order  of  decay  of  $(t),  to  the  right  hand  side. 


4.  foiasi liner  Problems 


We  investigate 

(4.1)  ey'  =  taA(z,t)v  +  taf(z,y,t,e) 

(4.2)  z'  =  tag(z,y,t,e) 


at  >  -1 
1  <  t  <  • 


«-3>  .*<*>( £!;.!)  -*«• 

(4.4)  (J)  e  C(  n,-] ) 


where  A(z,t)  is  an  n*n-matrix,  f  an  n-vector,  g  an  m-vector  and  we 
assume  that  the  Problem  (4.1),  (4.2)  is  quasilinear: 

(4.5)  57=°<e> 

for  t  €  f1,*»l,  e  6  [0feQ]  and  y,z  in  compact  sets.  We  get  immediately 

(4.6)  f (z,y,t,0)  =  f (z,0,t,0). 

We  now  assume  that  F(e)  is  a  k  *  (n+m)-matrix  (k  will  be  specified  later) 
and 


(4.7)  f,g  e  C2(Rm+n  X  (1,-1  x  [0,eQ))  nC^([1,-)) 

Now  we  proceed  as  Ringhofer  (1981)  did.  We  split  F(0)  into 

5 

(4.8)  F{ 0)  (^  )  «  Fz(°^2  +  Fy(0)5,,  5,  6  Rn,  ^  6  R1" 
and  we  assume  that  there  is  an  integer  r+  <  n  such  that 

(4.9)  F  (0)5,  =  F  (0)5*  +  F  (0)5~ 

I  7  .  7 _ 


where 

(4.10) 

and  the  k  x  r 


e  r  +, 


matrix  F  (0)(k  >  r 

y 


r 

e  R  ",  r_  -  n-r+ 


is  assumed)  has  maximal  rank 


r 
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Therefore,  there  is  a  k  *  k  matrix  Z  such  that 


(4.11) 


ZF  (0)  « 

v  ~ 

)r 

,  z  - 

y_ 

0 

r 

)k-r_ 

/ 

_  Z2. 
K~k 

holds  where  V  is  nonsingular. 

The  main  assumption  is  the  following.  The  reduced  problem 

(4.12)  z'  -  t°g(y,z,t,0),  1  <  t  <  • 

(4.13)  0  -  A(z,t)y  +  f(z,0,t,0),  1  <  t  <  ® 

(4.14)  Z2F(0)^z(lP  “  z20(°» 

(4.15)  (J)  e  C(H,**1) 

has  an  isolated  solution  (see  Keller  (1975))  y  *■  y,  z  *  z  and 


(4.16) 


A(z,t) 


holds  in  C  *  {(z,tlMz-S(t)l  <  <p,  t  6  [1,-1 ) ,  <P  >  0,  where  the 
eigenvalues  ^+(zft)  of  A+(z#t)  and  X^(zft)  of  A_(zft)  fulfill 

(4.17)  Re  X+(z,t)  >  c+  >  0  ,  (z,t)  e 

(4.18)  Re  Mz,t)  <  -c_  <  0,  (z,t)  e  C^. 

This  guarantees  that  (4.13)  can  be  solved  for  z  locally  around  z  and 

(4.19)  y  -  y(z,t)  -  -A_1(Eft)f (z,0,t,0 ) 
holds.  We  now  assume  that  the  matrix 
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(4.20) 
•1 


5  =  J2ur),s<“),“, 0)  +  J^U(-), y(-),-, 0)^f  <£<-),-). 

3  z 


(A  is  supposed  to  be  smooth)  fulfills 


(4.21) 


EJE 


-1 


+ 


J 


0 


where  the  eigenvalues  of  J+(J“)  have  positive  (negative)  real  parts. 

MW 

Therefore  we  assume  that  ZjFfO)  is  a  matrix  r_  x  (n+m)  matrix, 
r_ 

Z20(O)  e  R  and  k  =  r_  +  r_  such  that  (4.12),  (4.13),  (4.14),  (4.15)  is 
well  posed  with  respect  to  the  number  of  'finite*  boundary  conditions  (see 
Markowich  (1980a),  de  Hoog  and  Weiss  ( 1980a, b)).  Obviously  zm  **  z(®), 
y^,  «  y(“)  are  solutions  of 

(4.22)  (a)  0  -  g(yo.,zc,,-,0),  (b)  0  =  A(z0D,“)yn>  +  f(sl,,0,'M)) 

and  we  assume  that  z  ,y  are  isolated  and  that  f (y  ,z  ,t,0)  =  0, 

OO  OO  00  00 

9(y.»*.#t,0)  5  0,  t  >  5  >  1  holds.  Therefore  5  as  in  (4.20)  can  be 
calculated  a  priori  at  these  roots. 

Let  ♦( t,e )  denote  the  fundamental  matrix  of 

(4.23)  ev'  -  t°A(z(t),t)v,  *(1,6)  -  i. 

We  only  state  the  existence  result  since  the  proof  goes  along  the  lines  of  the 
proof  given  in  Ringhofer  (1981)  for  finite-interval  problems  using  the  linear 
theory  developed  in  chapters  2,3  of  this  paper. 

Theorem  4.1.  Let  F<6)  e  C([0,e  })  be  a  (r_+r_)  *  (n+m)  matrix.  Under  the 
given  assumption  the  problems  (4.1),  (4.2),  (4.3),  (4.4)  has  a  locally  unique 
solution  y,z  for  e  sufficiently  small  such  that 


yU,e) 


♦<t,e) 


C 


+  y(t)  +  0(e) 


z(t,e)  -  z(t)  +  0(6) 
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m 


H 


for  some  C  e  C  holds  uniformly  in  [I#*]* 
From  chapter  3  we  conclude  that 


l  +  lt^) 


01 


i  <  const.  exp(-£-^j^t  °+-1)),  c  >  0 


holds,  t-asymptotics  for  z(t),y(t)  can  be  obtained  from  Markowich  (1980a): 

(4.25)  z(t)  =  z(-)  +  E*<t,rjn*  +  oCM(t,r.°  i ,2) 


1  o  - 

C  +  o(i<Mt> 

“O' 

i~ 

Lr-j 

r_ 

for  Sec  where 

(4.26) 

holds.  From  (4.19)  we  get 

(4.27) 


♦  (  t )  =  exp(j~-(ta+1-1)) 


y(t)  *  y(«)  +  (i(»),-)E*(t)  x°ls  +  o(l*(t "|l2). 


The  approximating  'finite'  problems  are 

a  a 

(4.28)  *  ‘  ‘  ‘  '  ‘ 


£y;  31  t  A(zT,t)yT+t  f(zT,yTrt,e) 


1  <  t  <  T 


(4.29) 


■  t  9<yT«Vt,E) 


(4.30) 


ryT(1'EN 


(4.31) 


,yT(T'e> 

S(r-C,(ZT(T,C, )  *  l"’*' 


where  S(T,e)  is  a 


<W 


*  (n+m)  matrix  and  Y(T,e)  e  R 


r++  r+ 


We  choose 


(4.32) 


and 

(4.33) 


S  =  S(T,e ) 


0] 


0 


-[i  #o]&  (£<-),•) 

r+  3z 


Y(T,e)  =  s 


y(“) 
z(- ) 


Then  we  obtain 


(4.34) 


-  Tl 

z(T,e )' 


and  by  using  the  linear  stability  result  (2.54)  we  get  by  proceeding  as 
de  Hoog  and  Weiss  (1980a)  did 


(4.35) 


,  yT  -  y , 

■  LT  _  Ji 
ZT  z 


[1,T] 


<  — *  ■<;::!>  - *■ 


for  the  locally  unique  solution  yT,  z-  of  (4.28),  (4.29),  (4.30),  (4.31) 
such  that  (4.34)  constitutes  the  convergence  estimate. 

As  in  the  linear  case  this  asymptotic  boundary  condition  only  depends  on 
the  reduced  'infinite'  problem. 
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zero)  for  the  fast  and  slow  components  uniformly  on  the  whole  infinite  line. 
The  second  goal  of  the  paper  is  to  derive  convergence  estimates  for  the 
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20.  Abstract  (continued) 

solutions  of  'finite'  singular  perturbation  problems  obtained  by  cutting  the 
infinite  interval  at  a  finite  (far  out)  point  and  by  substituting  appropriate 
additional  boundary  conditions  at  the  far  end.  Using  a  suitable  choice  for 
these  boundary  conditions  the  order  of  convergence  is  shown  to  depend  only  on 
the  decay  property  of  the  infinite  solution. 


